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Using a recent, novel Hamiltonian formulation of the gravitational interaction of spinning binaries, 
we extend the Effective One Body (EOB) description of the dynamics of two spinning black holes to 
next-to-leading order (NLO) in the spin-orbit interaction. The spin-dependent EOB Hamiltonian is 
constructed from four main ingredients: (i) a transformation between the "effective" Hamiltonian 
and the "real" one, (ii) a generalized effective Hamilton- Jacobi equation involving higher powers of 
the momenta, (iii) a Kerr-type effective metric (with Pade-resummed coefficients) which depends 
on the choice of some basic "effective spin vector" S e ff, and which is deformed by comparable-mass 
effects, and (iv) an additional effective spin-orbit interaction term involving another spin vector er. 
As a first application of the new, NLO spin-dependent EOB Hamiltonian, we compute the binding 
energy of circular orbits (for parallel spins) as a function of the orbital frequency, and of the spin 
parameters. We also study the characteristics of the last stable circular orbit: binding energy, orbital 
' frequency, and the corresponding dimensionless spin parameter Slso = cJlso/(G(-Hlso/c 2 ) 2 ). We 

find that the inclusion of NLO spin-orbit terms has a significant "moderating" effect on the dynam- 
, 1 ^ 1 ical characteristics of the circular orbits for large and parallel spins. 

PACS numbers: 04.25.-g, 04.25.Nx 
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Coalescing black hole binaries are among the most promising sources for the currently operating ground-based 
network of interferometric detectors of gravitational waves. It is plausible that the first detections concern binary 
systems made of spinning black holes, because (as emphasized in [lj) the spin-orbit interaction can increase the 
binding energy of the last stable orbit, and thereby lead to larger gravitational wave emission. This makes it urgent to 
have template waveforms accurately describing the gravitational wave emission of spinning binary black holes. These 
waveforms will be functions of at least eight intrinsic real parameters: the two masses m\ , 7712 and the two spin vectors 
Si, S2. Due to the multi-dimensionality of the parameter space, it seems impossible for state-of-the-art numerical 
simulations to densely sample this parameter space. This gives a clear motivation for developing analytical methods 
for computing the needed, densely spaced, bank of accurate template waveforms. 

Among existing analytical methods for computing the motion and radiation of binary black hole systems, the most 
complete, and the most promising one, is the Effective One Body (EOB) approach [H, Q, H, 3 . This method was the 
first to provide estimates of the complete waveform (covering inspiral, plunge, merger, and ring-down) of a coales cing 
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INTRODUCTION 
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black hole binary, both for non-spinning systems [3( , and for spinning ones [5| . Several recent works (fj, [7|, la [9|, 
have shown that there was an excellent agreement 1 between the EOB waveforms (for non-spinning systems) and the 
results of recent numerical simulations (see ll| for references and a review of the recent breakthroughs in numerical 
relativity). In addition, the EOB method predicted, before the availability of reliable numerical relativity (NR) results, 
a value for the final spin parameter afi n of a coalescing black hole binary [U, Q which agrees within ~ 10% with the 
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1 For instance, Ref. [E| finds a maximal dephasing of ±0.005 gravitational wave cycles between EOB and numerical relativity waveforms 
describing 12 gravitational wave cycles corresponding to the end of the inspiral, the plunge, the merger and the beginning of the ringdown 
of an equal-mass coalescing binary black hole. 
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results of recent numerical simulations (see [ll[ for a review and references). Recently, it has been shown that the 
introduction of some refinements in the EOB approach, led to an EOB/NR agreement for Sfi n at the 2% level fl2| . 

In a previous paper [l[ the EOB method (originally developed for non-spinning systems) has been generalized to 
the case of spinning black holes. It was shown there that one could map the third post-Newtonian (3PN) orbital 
dynamics, together with the leading order (LO) spin-orbit and spin-spin dynamical effects of a binary system onto an 
"effective test particle" moving in a Kerr-type metric. In the present paper, we extend and refine the EOB description 
of spinning binaries by using a recently derived [l3| Hamiltonian description of the spin-orbit interaction valid at the 
next to leading order (NLO) in the PN expansion. (The NLO spin-orbit effects in the harmonic-gauge equations of 
motion were first obtained in [3 EH-) Let us recall that LO spin-orbit effects are proportional to G/c 2 , while NLO 
ones contain two sorts of contributions: oc G/c 4 and oc G 2 /c 4 . Regarding the spin-spin coupling terms, we shall use 
here only the LO results which arc made of two different contributions: the LO S1S2 terms (Tq| (which have been 
recently extended to NLO in [I3|), and the LO Sf and Sf terms. The latter are specific to Kerr black holes, being 
related to the quadrupole gravitational moment of a rotating black hole. 2 It was shown in [l| that the complete LO 
spin-spin terms (the sum of S1S2, Sf , and 5| terms) admitted a remarkable rewriting involving a particular linear 
combination So, defined below, of the two spin vectors. This fact, together with the more complicated structure of 
spin-orbit terms at the NLO, will lead us below to define a particular, improved EOB description of spinning binaries. 

The present paper consists of two parts: In the first part (Sections 2 and 3) we shall develop the formalism needed 
to finally define (in Section 4) our improved EOB description of spinning binaries. In the second part (Section 5), we 
shall consider one of the simplest "applications" of our EOB Hamiltonian: a discussion of the energetics of circular, 
equatorial orbits for systems with parallel spins. In this section, we shall make contact with previous related analytical 
investigations, notably [l5j . and prepare the ground for making contact with numerical data. 

A few words about our notation: We use the letters a, b = 1,2 as particle labels. Then, m a , x a = p a = (p a i), 
and S a = (S a i) denote, respectively, the mass, the position vector, the linear momentum vector, and the spin vector 
of the ath body; for a ^ b we also define r ao = x a — x , r a b = |r a i|, n ao = r ao /r aD , | • | stands here for the Euclidean 
length of a 3-vector. 



II. PN-EXPANDED HAMILTONIAN 



Our starting point is the PN-expanded (or "Taylor-expanded" ) two-body Hamiltonian H which can be decomposed 
as the sum of: (i) an orbital part H , (ii) a spin-orbit part _ff so (linear in the spins), and (iii) a spin-spin term H ss 
(quadratic in the spins), 

H(K a ,p a ,S a ) = H (x a ,p a ) + H so (x a ,p a ,S a ) 

+ ff ss (x a ,p a ,S Q ). (2.1) 

The orbital Hamiltonian H Q includes the rest-mass contribution and is explicitly known (in ADM-like coordinates) 
up to the 3PN order [3, [l9| . Its structure is 

H (x a ,p a ) = ^m a c 2 +iJ oN (x Q ,p a ) 

a 

+ \ ff lPN(x a , p a ) + -r -ffo2PN(x a ,Pa) 

+ ^5#o3PN(x o ,p o ) + Hg) • (2-2) 

The spin-orbit Hamiltonian H so can be written as 

) = ^n a (x D ,p D ).S Q , (2.3) 

a 

Here, the quantity Q a is the sum of a LO contribution (oc 1/c 2 ) and a NLO one (oc 1/c 4 ), 

fi o (x , P5 ) = ^°(X6, Pb) + "a L °(x D , Pi). (2.4) 



2 Note in passing that, if one wishes to describe the dynamics of, say, neutron-star binaries with the EOB formalism, one should add 
"correcting" S 2 and S 12 terms. 
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The 3-vectors and f2^ LO were explicitly computed in Ref. [l3j]. They are given, for the particle label a = 1, by 
fi i° = ^fT^ n i2xPi-2ni2xp 2 Y (2.5a) 




m\\ („ 15 



5— Ini2 x p! + I 6mi + ym 2 )ni 2 x p 2 



G i(_ 5m, 2 Pi _ 3(pi ■ p 2 ) 3p^ _ 3(ni2 ■ Pi)(ni 2 • p 2 ) _ 3(ni 2 ■ p 2 ) 2 
c4r i2 \ V Am\ 4mim 2 Am\ 2mim 2 



,' (pi-p 2 ) 3(ni2 -Pi)(n 12 -P2) \ v „ ,f 3(n 12 - pi) 2(n 12 ■ gg) , , , 
+ 1 ni2 x p 2 + — J Pi x p 2 ■ (2.5b) 



mim2 mii7i2 ) \ Amf mii7i2 

The expressions for and n% ho can be obtained from the above formulas by exchanging the particle labels 1 and 2. 

Let us now focus our attention on the dynamics of the relative motion of the two-body system in the center- of -mass 
frame, which is defined by the requirement Pi + p 2 = 0. It will be convenient in the following to work with suitably 
rescaled variables. We rescale the phase-space variables R = xi x 2 and P = pi = — p 2 of the relative motion as 
follows 

R Xi - x 2 P Pi P2 

p= — = — = , (2.6) 



GM GM ' ^ n fx [i 

where M = mi + m 2 and /i = mim 2 /M. Note that this change of variables corresponds to rcscaling the action by 
a factor l/{GMp). It is also convenient to rescale the original time variable T and any part of the Hamiltonian 
according to 

where _ff NR = H — Mc 2 denotes the "non relativistic" version of the Hamiltonian, i.e. the Hamiltonian without the 
rest-mass contribution. It has the structure H NR = -ip 2 — - + O (p-)- 

It will be convenient in the following to work with the following two basic combinations of the spin vectors: 

S = Si + S 2 = rriicai + 77i2ca 2 , (2.8a) 
S* = — -Si H ^-S 2 = m 2 cai + mica 2 , (2.8b) 

771 1 777 2 

where we have introduced (as is usually done in the general relativistic literature) the Kerr parameters 3 of the 
individual black holes, ai = Si/(mic) and a 2 = S 2 /(m 2 c). Note that, in the "spinning test mass limit" where, say, 
m 2 — ► and S2 — > 0, while keeping a 2 = S , 2 /(t71 2 c) fixed, we have a "background mass" M ~ mi, a "background 
spin" Sbckgd = Mcabckgd — Si = micai, a "test mass" /1 ~ 77i 2 , and a "test spin" Stest = S 2 = 77i 2 ca 2 ~ /icatest 
[with a tos t = Stest /(m c )]- Then, in this limit the combination S ~ Si = m\c&.\ ~ Mc abekgd = Sbckgd measures 
the background spin, while the other combination, S* ~ mica2 ~ Mca tes t = AfStest/A* measures the (specific) test 
spin atest = Stest /(/^c). The quantities S and S* are the two simplest symmetric (under the permutation 1 <-> 2) 
combinations of the two spin vectors which have these properties. 

In view of the rescaling of the action by a factor l/(GM/i), corresponding to the rescaled phase-space variables 
above, it will be natural to work with correspondingly rescaled spin variables 4 



3 Note that we use here the usual definition where the Kerr parameter a = S/ (Mc) has the dimension of length. We denote the associated 
dimensionless rotational parameter with an overhat: a = ac 2 /(GM) = c S/ (GM 2 ). 

4 We recall that (orbital and spin) angular momenta have the same dimension as the action. 
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for any label X (X = 1,2, *,•••)• 

Making use of the definitions (|2.6p ~ (|2.9[) one easily gets from Eqs. (|2.3|) - (|2.5p the center-of-mass spin-orbit Hamil- 
tonian (divided by /z) expressed in terms of the rescaled variables: 

A ; a g so (r,p,S,S*) 
H so (r, p, S, S ) = 

A* 

= Ii?£°o(r,p,S,S*) 

+ i^ (r,p,S,S*) + o(i], (2.10) 



where (here n = r/|r|) E 



iJ£°o(r,p,S,S*) = ^l2{S,n,p) + |(S*,n,p) I, (2.11a) 



i^° LO (' , P, S, S* ) = ^ i - (6 + 2i/) (5, n, p) - (5 + 2v) (S* , n, p) 



£l (y^P 2 + ^( n 'P) 2 )(^^p) 



| + 2^ p 2 + \v (n • p) 2 ) (5*, n,p) 1 , (2.11b) 



with v = n/M ranging from (test-body limit) to 1/4 (equal-mass case). 
Note that the structure of the rescaled spin-orbit Hamiltonian is 



tf so (r,p,S,S*) = JL;(g*™(S,n,p) +gj? u (S*,n,p)). (2.12) 
This corresponds to an unrescaled spin-orbit Hamiltonian of the form 

G L / ADM adm c *\ ro ti\ 

Hso = ^2^3 ' {9s s +.9s* s J. ( 2 - 13 ) 

where R = GMr is the unrescaled relative distance (in ADM coordinates), LeRxP = GM/ir x p the relative 
orbital angular momentum, and where we have introduced two dimensionless coefficients which might be called the 
"gyro-gravitomagnetic ratios" , because they parametrize the coupling between the spin vectors and the "apparent" 
gravitomagnetic field 

GM RxP 

V X V — OC 

c 2 R R 3 

seen in the rest-frame of a moving particle (see, e.g., Refs. 0, Hl| for a discussion of the expression of the "grav- 
itomagnetic field" in the rest-frame of a moving body). The explicit expressions of these two gyro-gravitomagnetic 
ratios are 

ff ADM = 2 + -L + \v (n • p) 2 - (6 + 2v) I) , (2.14a) 

^((-^ + 2,)p 2 + ^(n.p) 2 -( 5 + 2,)i). (2.14b) 



fJ's 



ADM _ ^ 



We introduce the following notation for the Euclidean mixed product of 3-vectors: [Vi, Vi, V3) = Vi ■ (V2 x V3) = Sij^V^V^V^ 
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In the following we shall introduce two related "effective" "gyro-gravitomagnetic ratios" , that enter the effective 
EOB Hamiltonian (in effective coordinates). The label "ADM" on the gyro-gravitomagnetic ratios (|2.14p is a reminder 
of the fact that the NLO value of these ratios depend on the precise definition of the radial distance R (which is 
coordinate dependent). Let us, however, briefly discuss the origin of the (coordinate-independent) LO values of these 
ratios, namely 

«§° = 2, 9^ = \=2-\. (2.15) 

Here the basic ratio 2 which enters both and g^P comes from the leading interaction, predicted by the Kerr 
metric, between the orbital angular momentum of a test particle and the background spin. See Eq. (|4.17|) below. As 
for the — i "correction" in the coupling of the "test mass" spin combination S* it can be seen (e.g. from Eq. (3.6b) 
of (22J) to come from the famous | factor in the Thomas precession (which is a universal, special relativistic effect, 
separate from the effects which are specific to the gravitational interaction, see Eqs. (3.2) and (3.3) in (22I]). 

To complete this Section, let us recall the remarkable form [found in Ref. see Eq. (2.54) there] of the leading- 
order spin-spin Hamiltonian H ss (including Sf, Sf as well as S1S2 terms). The unrescaled form of the spin-spin 
Hamiltonian reads 

H aB (R i S ) = ~S l S 3 Q a ij ^, (2.16) 



while its rescaled version reads 



H 

1^5* & d .l ^ 2 3(n-S )^Sg 
"2 c 2 lJ r~ 2 c 2 r3 ' [ ' 

The remarkable fact about this result is that it is entirely expressible in terms of the specific combination of spins 
So = GM/iSq defined as: 

So ee S + S* = (l + ^) Sl + f 1 + (2.18) 

We shall come back below to the remarkable properties of the combination So, which will play a central role in our 
EOB construction. 



III. EFFECTIVE HAMILTONIAN AND "EFFECTIVE GYRO-GRAVITOMAGNETIC" RATIOS 

We have obtained in the previous Section the expression of the full center-of-mass- frame Hamiltonian (|2.1[) , in PN- 
expanded form. In order to transform this Hamiltonian into a format which can be resummed in a manner compatible 
with previous work on the EOB formalism, we need to perform two operations on the Hamiltonian (|2.ip . First, we 
need to transform the phase-space coordinates (x a ,p a ,S Q ) by a canonical transformation compatible with the one 
used in previous EOB work. Second, we need to compute the effective Hamiltonian corresponding to the (canonically 
transformed) real Hamiltonian (|2.1[) . 

We start by performing the purely orbital canonical transformation which was found to be needed in Refs. 0, 3 to 
go from the ADM coordinates (used in the PN-expanded dynamics) to the coordinates used in the EOB dynamics. 
This orbital canonical transformation is (implicitly) given by 

a/W + ^A K= Pi -^#i (3.1) 
op\ ax 1 

Here the orbital generating function G (q,p') has been derived to 2PN accuracy in and to 3PN accuracy in pj. In 
the present paper, as we are only concerned with the additional spin-orbit terms, treated to 1PN fractional accuracy, 
it is enough to work with the lPN-accurate generating function G (x,p'). In terms of the rescaled variables, the 
rescaled lPN-accurate orbital generating function reads 

- G (r,p) 
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= ? (r.p)^- 5 ^+^l + r j-J. (3.2) 

This transformation changes the phase-space variables from (r, p, S, S*) to (r', p', S, S*). At the linear order in the 
transformation (which will be enough for our purpose), the effect of the transformation on any of the phase-space 
variable, say y, is y' = y + {y, G }, where {•, •} denotes the Poisson bracket. As G is independent of time, it leaves the 
Hamiltonian numerically invariant: H'(y') = H{y). This means that it changes the functional form of the Hamiltonian 
according to H'{y') = H(y' — {y, G }) = H(y') — {H, G }. Note the appearance of the opposite sign in front of the 
Poisson bracket, with respect to the effect of the generating function on the phase-space variables. 

As G is of order 1/c 2 , its explicit effect on the two separate terms, Hf^ Q and H^ LO , in the PN expansion of the 
spin-orbit Hamiltonian is given by: 

^ s o°(r',p',S,S*) =#fo(r',p',S,S*), (3.3a) 

^ s L °o(r', p', S, S*) = H s N l (v', p', S, S*) 

-{ff£°o,G }(r',p',S,S*). (3.3b) 

It will be convenient in the following to further transform the phase-space variables by performing a secondary purely 
spin-dependent canonical transformation, affecting only the NLO spin-orbit terms. The associated new generating 
function, G s (r, p, S, S*) (assumed to be proportional to the spins and of order 1/c 4 ) will change the variables (y') = 
(r', p', S, S*) into (y") = (r", p", S", S"*) according to the general rule 6 y" = y' + {y', G s }. For the same reason as 
above, the (first-order) effect of G s on the functional form of the Hamiltonian will involve a Poisson bracket with the 
opposite sign: H"{y") = H{y") - {H,G S }. 

We shall consider a generating function whose unrescaled form reads 



G S (R, P, S, S*) = A ig (R • P)(R x P) • (a{v) S + b{v) S*) , (3.4) 
while its rescaled form reads 

G s (r,p,S,S*) 



_ G S (R,P,S,S*) 




a(u)(S,n,p)+b(v)(S*,n,p)]. (3.5) 



Here a(y) and b(v) are two arbitrary, ^-dependent dimcnsionless coefficients. 7 Similarly to the result above, the 
explicit effect of this new canonical transformation on the two separate terms, ■Hlo an< ^ -^nlc m ^ ne expansion 
of the spin-orbit Hamiltonian reads: 

fff&V'.p", S", S"*) - ^ s °(r",p", S", S"*), (3.6a) 

TTll SO /_// // Bl/ 0"*^ TJ ,so (r- 11 n" «" 

-"NLO^ >P ' 3 ' 3 I — -"NLOl r iP i s i s ) 

-{ff oN ,G s }(r",p",S",S"*), (3.6b) 

where _ff N is the Newtonian orbital Hamiltonian. In the following, we shall, for simplicity of notation, omit the 
double primes on the new phase-space variables (and on the corresponding Hamiltonian). 

The second operation we need to do is to connect the "real" Hamiltonian H to the "effective" one H e g , which is more 
closely linked to the description of the EOB quasi-geodesic dynamics. The relation between the two Hamiltonians is 
quite simple @, Q : 



H c s _ H 2 - mfc 4 - m\c A 
fie 2 2miTO2C 4 



, , ( 3 - 7 ) 



6 Note that while Go did not affect the spin variables, the spin-dependent generating function G s will now affect them. 

7 The coefficients a(u) and b(u) can be thought of as being two "gauge" parameters, related to the arbitrariness in choosing a spin- 
supplementary condition, and in defining a local frame to measure the spin vectors. 
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where we recall that the real Hamiltonian H contains the rest-mass contribution Mc 2 = (mi + m-^c 1 . Let us also 
note that Eq. f|3.7|) is equivalent to 



,2 =1 + -72T + ^^7274^' ( 3 - 8 ) 



H cS H NR 1 <H NR ) 2 
— ^- = H h -1/- — 5- 

/xc 2 /ic 2 2 ^ 2 c 

where _ff NR denotes the "non relativistic" part of the total Hamiltonian H, i.e., i? NR = H — Afc 2 , or more explicitly 



rjNR /rr , ^olPN . -ffo2PN . #o3PN' 
— I -HoN H n 1 -. 1 ; — 



zjso t/so 
^LO 1 "NLO 



(3.9) 



By expanding (in powers of 1/c 2 and in powers of the spins) the exact effective Hamiltonian (|3.7[) . one easily finds 
that the "spin-orbit part" of the effective Hamiltonian H e ff (i.e. the part which is linear-in-spin) differs from the 
corresponding part H so in the "real" Hamiltonian by a factor ~ 1 + vH Nn /c 2 ~ 1 + vH ^/c 2 , so that we get, for the 
explicit PN expansion of H^, 

U so 1 1 / \ 

-~ = Tj-^lo + 774 (-^nlo + vH oN Hl° ) . (3.10) 

Combining this result with the effect of the two generating functions discussed above (and omitting, as we already 
said, the double primes on the new phase-space variables (r", p", S", S"*)), we get the transformed spin-orbit part of 
the effective Hamiltonian in the form 



rrso 



which corresponds to the following unrescaled form (with L = R x P 



(nxp).( g fS+. 9 ^S*), (3.11) 



^ ff = §^-fe ff S + ^S*). (3.12) 



c 



Here the two "effective gyro-gravitomagnetic" ratios g^ 1 and g^l differ from the "ADM" ones introduced above by 
three effects: (i) a factor ~ 1 + vH NR /c 2 ~ 1 + vH Q ^/c 2 due to the transformation from H to H e ff, (ii) the effect of 
the orbital generating function G a going from ADM to EOB coordinates, and (iii) the effect of the spin-dependent 
generating function G s , which involves the gauge parameters a(v) and b(y). Their explicit expressions are then found 
to read 



gf = 2 + 1 + a(u)) p 2 + 3a(v)) (n • p) 2 - i (v + a{v))\ (3.13a) 

gf, = I + I (( - \ + \v + b{v)) P 2 - (H v + 36(1/)) (n • P ) 2 - I (I + \ v + 6(,))) . (3.13b) 

The choice of the two "gauge" parameters a(v) and b(v) is arbitrary, and physical results should not depend on 
them. 8 This would be the case if we were dealing with the exact Hamiltonian. However, as we work only with an 
approximation to the exact Hamiltonian, there will remain some (weak) dependence of our results on the choice of 
a(v) and b(v). We can use this dependence to try to simplify, and/or to render more accurate, the spin-orbit effects 
implied by the above expressions. In particular, we shall focus in this paper on a special simplifying choice of these 
gauge parameters: namely, the values 

a{v) = -K, &(,/) = jj-i„, (3.14) 

which suppress the dependence of the effective gyro-gravitomagnetic ratios on p 2 . With this particular choice, the 
explicit expressions of these ratios become 



Note in particular that the gyro-gravitomagnetic ratios do not depend on a(y) and b(y) when considering circular orbits, i.e. when 
p 2 = 1/r and (n ■ p) = 0. 
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gf = 2 + ±(~v(n.p) 2 --,<-). C-!.:i5a.! 



5 l\ 

rf-i^(-(T + 5^-(S^0F> ^ 

IV. SPIN-DEPENDENT EFFECTIVE-ONE-BODY HAMILTONIAN 

Up to now we only considered PN-expanded results. In this Section, we shall generalize the approach of [l[ in 
incorporating, in a resummed way, the spin-dependent effects within the EOB approach. Let us first recall that the 
approach of [l| consists in combining three different ingredients: 

• a generalized Hamilton- Jacobi equation involving higher powers of the momenta (as is necessary at the 3PN 
accuracy [4]); 

a/3 



a z^-deformed Kerr- type metric g c g , which depends on the choice of some basic "effective spin vector" S^g', 



• the possible consideration of an additional spin-orbit interaction term A7J so (r, p, So, cr) in the effective Hamil- 
tonian, whose aim is to complete the spin-dependent interaction incorporated in the definition of the Hamilton- 
Jacobi equation based on a certain choice of "effective spin vector" S\ s . 

At the LO in spin-orbit and spin-spin interactions, Ref. [l[ showed that one had the choice between two possibilities: 

(i) use as effective spin vector the combination S + |S* which correctly describes the LO spin-orbit effects, but only 
approximately describes the LO spin-spin effects; 9 or 

(ii) use as effective spin vector the combination 

S .S + S* = fl + -H + fl + -K (4-1) 



which correctly describes the full LO spin-spin interaction (see (|2.I7p above), and complete the description of the LO 
spin-orbit effects by adding a term AH ao (r, p, So, <x) involving a suitably defined spin combination cr. (At LO, Ref. 
Q defined er LO = -|S*.) 

Intuitively speaking, the second possibility consists in considering that the "effective particle" is endowed not only 
with a mass [i, but also with a "spin" proportional to cr, so that it interacts with the "effective background spacetime" 
both via a geodesic- type interaction (described by the generalized Hamilton- Jacobi equation), and via an additional 
spin-dependent interaction proportional to its spin oc cr. 

At the present, NLO approximation, where it is crucial to accurately describe the spin-orbit interaction, as well as, 
by consistency, the LO spin-spin ones, we have chosen to follow the second possibility, which offers more flexibility, 
and which looks natural in view of the remarkably simple LO result (|2.17|) for the spin-spin interaction (see, however, 
the suggestion at the end of the concluding Section 6). 

Therefore we shall successively introduce the ingredients needed to define 

• the Hamilton- Jacobi equation (describing the basic "geodesic- type" part of the effective Hamiltonian); 

• the effective, z/-deformed Kerr- type metric g"£ ; 

• the "effective spin vector" 5* ff entering the previous Kerr-type metric; 

• the additional spin-orbit interaction A_ff so (r, p, So, cr) involving a new, specific NLO spin combination cr. 
The modified Hamilton-Jacobi equation [i| is of the form 

9 a l P a Pp + Qi(Pi) = -M 2 c 2 , (4.2) 



9 One can then correct for the missing terms by adding an explicit supplementary term in the Hamiltonian, quadratic in the spins. 
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where Qi(Pi) is a quartic-in-momenta term (which only depends on the space momentum components Pi). For 
circular orbits Qi{P{) will be zero (see 0, so that we will not need its explicit expression in the present paper. 

The role of the Hamilton- Jacobi equation above is to allow one to compute the main part (modulo the additional 
spin-orbit interaction added later) of the effective Hamiltonian H"^ ln — E e g = —Pqc by solving (|4.2p with respect to 
P . The result can be written as 



jjniain 
^cff 



E, 



off 



PPiC- 



/I 2 C 2 



7 ij PiPj+Q4(Pi), 



where we have introduced the auxiliary notation 



n = (-a 00 )- 1 / 2 = V-J' = n v 9cS 9cS 

a — \ 9cs) > P — oo ' 7 — 9 c s 



9° C ° S 



(4.3) 



(4.4) 



The next crucial ingredient consists in defining the (spin-dependent) effective metric entering the Hamilton- Jacobi 
equation, and thereby the effective Hamiltonian (|4.3[) . We shall follow here Ref. [I[ in employing an effective co-metric 
of the form (here Pt = cPo) 



P lb + asm 2 9—^ 2 



1 



MR) 



(R 2 + a 2 )^+aP^) \, 



where the functions A t and are defined as 

A t (R) = R 2 P£ 



A(R) 



R 2 



A R (R) = A t (R)D- 1 (R), 



(4.5) 

(4.6a) 
(4.6b) 



and where the Kerr-like parameter a is defined as a = S c ft / (Mc) , where S e e denotes the modulus of the "effective 
spin vector" S** ff entering the definition of the Kerr-like metric above. We shall come back below to the choice of this 
vector Sl s (which is one of the ingredients in the definition of a spin-dependent EOB formalism) . In Eq. (|4.6a|) P^ 
denotes the operation of taking the (n,m)-Pade approximant, 10 and the PN expansions of the metric coefficients A 
and D~ x equal (here u = GM/{Rc 2 )) 



/94 41 \ 
A(u) = l-2u+ 2//ti 3 + (y - 32 J 

D _1 (u) = I + &vu 2 + 2(26 - 2,v)vu z 



1 \ -4 
7T \VU , 



(4.7a) 
(4.7b) 



For pedagogical clarity, we have given above the expression of the effective EOB metric in a Boyer-Lindquist-type 
coordinate system aligned with the instantaneous direction of the (time-dependent) effective spin vector S* s . This 
expression will suffice in the present paper where we will only consider situations where the spin vectors are aligned 
with the orbital angular momentum, so that they are fixed in space. As emphasized in when applying the EOB 
formalism to more general situations (non aligned spins) one must rewrite the effective co-metric components in a 
"fixed" Cartesian-like coordinate system. This is done by introducing 



rt^x'/R, s l = ^, cos6 = nV, 

Jeff 

p = \J R 2 + a 2 cos 2 6, 



(4.8) 



10 Let us recall that the (n, m)-Pade approximant of a function cq + C1U + C2U 2 + ■ ■ ■ + c n _|_ m M"+ m is equal to N n (u) / D m (u) , where N„(u) 
and D m (u) are polynomials in u of degrees n and m, respectively. 
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and rewriting the co-metric components as 



a _ j* (4 , a) 



° 2 (sxR)'( S xR)J. (4.9c) 
At (it) 



Making use of Eqs. (|4.9[) one computes 



A*(i?) 



(i? 2 + a 2 ) 2 - a 2 A t (R) sin 2 9 
a(R 2 + a 2 - A t (R)) 



(4.10a) 



{R 2 + a 2 ) 2 - a 2 A t (R) sin 2 i 



(s x R)\ (4.10b) 



8* ft> 

7 y =^ + ^-- (4-lOc) 

Replacing the latter expressions in the general form of the effective energy (|4.3p yields the most general form of the 
main part of the effective Hamiltonian i/^ am (x, P, S a ). 

The definition of H"^ ln (x, P, S a ) crucially depends on the choice of effective Kerr-type spin vector. In order to 
automatically incorporate, in a correct manner, the LO spin-spin terms, we shall use here 

Mca= S cff ee S = S + S* = + — )Si + (l + — W (4.11) 

V nil ' ^ m 2 ' 

Note that, besides its usefulness in treating spin-spin effects, this definition has several nice features. For example, if 
we introduce the Kerr parameters of the individual black holes, ai = Si /(Mc), & 2 = S 2 /(Mc), the Kerr parameter 
ao = So /(Mc) (where we naturally take mo = M = mi + m-i) associated to the spin combination (|4.1j) is simply 

a =ai+a 2 . (4.12) 

Let us also note that the corresponding dimensionless spin parameters (with, again, m = M = mi + m%) 

^ = ^' * = 0,1,2, (4.13) 

satisfy 

a = Xiki +X 2 & 2 , (4.14) 

where Xi = mi/M and X 2 = m 2 /M are the two dimensionless mass ratios (with Xi + X 2 = 1 and XiX 2 = v). 
This last result shows that, in a-space, the "point" ao is on the straight-line segment joining the two "points" ai and 
a 2 . The individual Kerr bounds tell us that each point ai and a 2 is contained within the unit Euclidean sphere. By 
convexity of the unit ball, we conclude that the "effective" dimensionles spin parameter a will also automatically 
satisfy the Kerr bound |ao| < 1. This is a nice consistency feature of the definition of the associated Kerr-type metric. 

It remains to define the additional "test-spin" vector <r, and the associated additional effective spin-orbit interaction 
term. Following the logic of [3 (and generalizing the LO results given in Eqs. (2.56)-(2.58) there), these quantities 
are defined by 

° = \gfs + \gf*s* - Scff 

= l( g f-2)S + ±(gf.-2)S*, (4.15) 

and 
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Aii so (x, P, 5s 0) oj = ttp; 575 9 . . 2 77 , (4.1b) 

(R 2 + al) 2 - alAt(R) sin 2 6 M 

where a = So/(Mc) and cos#o = i*<So/|So|- The justification for these definitions is that the "main" Hamilton- Jacobi 
part of the effective Hamiltonian contains, as spin-orbit (i.e. linear-in-spin) part, the following term 



linear-in-spin 

R 2 + a 2 - A t (R) 
(R 2 +a 2 ) 2 -a 2 A t (R)sin 2 fl, 

2GM 



= cP z tLJL *w - (ao x R) 



linear- in- spin 



cR 3 
c 2 R 3 



Pi(a x R)* + (NNLO corrections) 

S + (NNLO corrections), (4.17) 



where the factor 2GM comes from the second term in the PN expansion of A t (R) — R 2 — 2GMR/c 2 + 
2v{GM) 3 / (Rc 6 ) + (quadratic-in-spin terms). Note that the absence of c -4 correction in the effective metric function 
A(R) means that the leading term oc 2GM in the spin-orbit part of H mmn is valid both to LO and to NLO, i.e., up 
to "next to next to leading order" (NNLO). 

When comparing this result to the NLO result (|3.12p . we see that the "main" part of the effective Hamiltonian 
contains a spin-orbit piece which is equivalent to having effective gyro-gravitomagnetic ratios equal to g™ aincff = 2 
and <7gl amcff = 2, instead of the correct values derived above. One then easily checks that the definition above of a 
and of the associated supplementary spin-orbit interaction A7J so (x, P, So, <r) has the effect of including the full result 
for the NLO spin-orbit interaction. It is also to be noted that the additional spin-orbit interaction AH so goes to zero 
proportionally to v in the test mass limit mi — > because, on the one hand, g c <? — 2 is proportional to v (if a{y) is), 
and, on the other hand, though g^l — 2 does not tend to zero with v, the second spin combination S* does tend to 
zero proportionally to v [see Eqs. (|5.2[) below]. 

Summarizing: we propose to define a total effective spin-dependent Hamiltonian of the form 

ff eff (x, P, Si, S 2 ) = ff e T n (x, P, S ) + Aff so (x, P, So, tr), (4.18) 

where iJ^ am (x, P, So) is given by the right-hand side of Eq. (|4.3p computed for the effective spin variable equal to S 
[defined in Eq. (|4.ip ] and where AiJ so (x, P, So, <x) is the additional spin-orbit interaction term defined above [with 
ao = S /(Mc)]. 

Finally, the real EOB-improved Hamiltonian (by contrast to the "effective" one) is defined by solving Eq. (|3.7[) with 
respect to H Tcai = H NR + Mc 2 : 



flrcal = Mc 2 Jl + 2u(^L -l), (4.19) 



where H e g is given in Eq. (|4.18[) . 



V. DYNAMICS OF CIRCULAR ORBITS 



In this Section we shall apply the construction of the NLO spin-dependent EOB Hamiltonian to the study of the 
dynamics of circular orbits of binary black hole systems. 

Besides the dimensionless spin parameters ai and a 2 already introduced above, it is convenient to introduce the 
dimensionless spin variables corresponding to the basic spin combinations S and S*, namely 

„ a „ a* 

(5.1) 



GM 2 ' GM 2 ■ 

Let us note in passing the various links between the dimensionless spin parameters that one can define [including 
a = cS a /(GM 2 ) already introduced above], 



k = Xlki + Xlk 2l k* = vk 1 + vk 2 , 



(5.2a) 
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a = a + a* = + X 2 a 2 . (5.2b) 

Here as above we use the mass ratios X\ = mi/M, X 2 = m 2 /M such that X\ + X 2 = 1 and X\X 2 = v. Let 
us note that for equal-mass binaries (mi = vn 2l X\ = X 2 = h), with arbitrary (possibly unequal) spins, one has 
a = a* = i(ai + a 2 ) = ^a . Note also that, in the test-mass limit, say mi 3> m 2 so that X% — ► 1 and X 2 — > 0, one 
has 

a = ao = a 1; a* = 0. (5.3) 
In the general case where the spin vectors are not aligned with the (rescaled) orbital angular momentum vector 11 

£, 

£ = r n x p, (5.4) 

there exist no circular orbits. However, there exist (at least to a good approximation) some "spherical orbits", i.e. 
orbits that keep a constant value of the modulus of the radius vector r, though they do not stay within one fixed 
plane. As discussed in [l[ one can analytically study these spherical orbits within the EOB approach, and discuss, in 
particular, the characteristics of the last stable spherical orbit. 

For simplicity, we shall restrict ourselves here to the situation where both individual spins are parallel (or antiparal- 
lel) to the orbital angular momentum vector £. In that case, we can consistently set everywhere the radial momentum 
to zero, p r = n • p = 0, and express the (real) EOB Hamiltonian as a function of r, £ = p v (using p 2 = £ 2 /r 2 , where 
I = |£|), and of the two scalars a, a* measuring the projections of our basic spin combinations on the direction of the 
orbital angular momentum £. They are such that 

a-£ = a£, a*-£ = a*£. (5.5) 

The scalars a and a* can be either positive or negative, depending on whether, say, a is parallel or antiparallcl to £. 
The sequence of circular (equatorial) orbits is then determined by the constraint 

dH rced (r,e,a,a*) 
dr 

Then, the angular velocity along each circular orbit is given by 

1 dH Ieal (r,£,a,a*) 

~ gmji m ■ (5J) 

As mentioned above, we have chosen the special values a{v) = — b(v) = | — ~v of the two gauge parameters, to 
simplify the expression of the Hamiltonian. 

In Figs. 1-4 we explore several aspects of the dynamics of circular orbits, using as basic diagnostic the relation 
between the energy and the angular velocity along the sequence of circular orbits ( "binding energy curve" ) . More 
precisely, we plot the dimensionless "non relativistic" energy 

as a function of the dimensionless angular velocity: 

Cl = ^-fl. (5.9) 

c 6 

For simplicity, we shall restrict most of our studies to symmetric binary systems, i.e. systems with mi = m 2 and 
ax — a 2 . For such systems the dimensionless effective spin parameter is eto = a\ = a 2 . The information contained in 
these figures deals with the following aspects of the description of the dynamics: 

• As a warm up, and a reminder, Fig. 1 considers the case of non-spinning binaries (i.e. So = 0). This figure 
contrasts the behaviour of the successive PN versions of the EOB dynamics, with that of the successive PN 
versions of the non-resummed, "Taylor-expanded" Hamiltonian. The numbers 1,2,3 refer to 1PN, 2PN, and 



In the following, we switch again to the use of scaled variables: r = R/(GM), £ = L/(GM/i), and p = P/fi. 
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FIG. 1: Binding energy curves for circular orbits of symmetric non-spinning binaries (mi = m,2 and ai = a2 = 0): dimensionless 
"non relativistic" energy e versus dimensionless angular frequency Cl. The notation E(n, *) means computation of the energy 
using the EOB-improved real Hamiltonian (|4.19|l with the nPN-accurate metric function At (R); the function At(i?) was 
computed by means of Eq. (|4,6a|l using the (1, n) Pade approximant at the nPN order. Here n = 1, 2, 3, 4, where n = 4 refers 
to the "4PN" case where a term -\-a$vu 5 is added to the function A(u). For the curves labelled by T(n, *) the computation was 
done with the direct PN-expanded (ADM-coordinates) orbital Hamiltonian (|2.2[l with the terms up to the nPN order included. 



3PN, while the letter "E" refers to "EOB" and the letter "T" refers to "Taylor" . For instance, E(3, *) refers 
to the e(f2) binding energy curve computed with the 3PN-accurate EOB Hamiltonian. [The star in E(3, *) 
replaces the label we shall use below to distinguish LO versus NLO treatment of spin-orbit effects. In the 
present non-spinning case we are insensitive to this distinction.] To be precise, the notation E(n, *) refers to 
a computation of the circular orbits using the do — > limit 12 of the EOB-improved real Hamiltonian (|4.19p 
with the nPN-accurate metric function A t (R); where A t (i?) was computed by means of Eq. (|4.6b[) using the 
following Pade approximants: (1,1) at the 1PN order, (1,2) at the 2PN order, and (1,3) at the 3PN order. As 
for the Taylor-based approximants to the binding energy curve, T(n, *), they were computed by using as basic 
Hamiltonian (to define the dynamics) the nPN-accurate Taylor-expanded Hamiltonian, in ADM coordinates, 
(12. 2|) , without doing any later PN re-expansion. 13 

It is interesting to note that the successive PN-approximated EOB binding energy curves are stacked in a 
monotonically decreasing fashion, when increasing the PN accuracy, and all admit a minimum at some value of 
the orbital frequency. This minimum corresponds to the last stable circular orbit (see below). The monotonic 
stacking of the EOB energy curves therefore implies that a higher PN accuracy predicts circular orbits which 
are more bound, and can reach higher orbital frequencies. Let us note in this respect that recent comparisons 
between EOB and numerical relativity data have found the need to add a positive 4PN additional term -\-a5vu 5 
in the basic EOB radial potential A(u) of Eq. (|4.7ap above, with 05 somewhere between +10 and +80 0, H, H, El . 
Though we do not know yet what is the "real" value of the 4PN coefficient 05 we have included in Fig. 1 two 
illustrative 14 values of this "4PN" orbital parameter, namely 05 = +25 and 05 = +60. Note that the effect of 
such positive values of 05 is to push the last few circular orbits towards more bound, higher orbital frequency 
orbits. This effect will compound itself with the effects of spin explored below, and should be kept in mind when 
looking at our other plots. 

By contrast with the "tame" and monotonic behaviour of successive EOB approximants, we see on Fig. 1 that 



Note that the So — > limit of the Pade resummation of some ao-dependent metric coefficient is not necessarily the same as the Pade 
approximant one might normally consider in the non-spinning case. 

As is well-known there are always many non-equivalent ways of defining any "nPN" result, depending of where, and how, in the 
calculation one is replacing a function by a PN-expanded polynomial. For instance, one could PN re-expand the function giving the 
energy e in terms of the orbital frequency Q, or the function giving e in terms of the orbital angular momentum L (see Ref. y for 
the computation of several such functions in the non-spinning case). However, we are ultimately interested (for gravitational- wave 
purposes) in defining a complete dynamics for coalescing spinning binaries. Therefore, we focus here on the results predicted by 
Hamiltonian functions H(x,p, • ■ ■ ). 

These two values of the 4PN parameter 05 were found in Refs. 0,13 t° f> e representative of the values of 05 that improve the agreement 
between EOB waveforms and numerical relativity ones. 
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FIG. 2: Binding energy curves for circular orbits of symmetric parallely spinning binaries (mi = mi and ai = ai oc r x p): 
dimensionless energy e versus dimensionless angular frequency Cl along circular orbits for various values of the dimensionless 
effective spin parameter So = cSo/{GM 2 ) = Si = 02 within the effective-one-body approach. The label E(3, 1) means that 
we use the EOB Hamiltonian with 3PN-accurate orbital effects and NLO spin-orbit coupling, i.e. Eq. (|4.15[1 was used with the 
NLO gyro-gravitomagnetic ratios g^ and g$l, Eqs. (|3.13fl . 

the successive Taylor-Hamiltonian approximants T(n, *) have a more erratic behaviour. Note in particular, that 
the 3PN-accurate energy curve does not admit any minimum as the orbital frequency increases (in other words, 
there is no "last" stable circular orbit). In view of this bad behaviour of the 3PN-accurate orbital Taylor- 
Hamiltonian, we shall not consider anymore in the following figures the predictions coming from such Taylor 
Hamiltonians. 15 

• In Fig. 2 we study the effect of changing the amount of spin on the black holes of our binary system. We use 
here our new, NLO spin-orbit EOB Hamiltonian, as indicated by the notation E(3, 1), where the first label, 3, 
refers to the 3PN accuracy, and the second label, 1, to the 1PN fractional accuracy of the spin-orbit terms (i.e., 
the NLO accuracy). Note that the EOB binding energy curves are stacked in a monotonically decreasing way 
as the dimensionless effective spin &q increases from ao = — 1 (maximal spins antiparallel to the orbital angular 
momentum) to ao = +1 (maximal spins parallel to the orbital angular momentum). Note also that this curve 
confirms the finding of |l| that parallel spins lead to the possibility of closer and more bound circular orbits. 

• Fig. 3 contrasts the effect of using the NLO spin-orbit interaction instead of the LO one in the EOB Hamiltonian. 
We use the full 3PN accuracy, and include the LO spin-spin interaction. E(3, 0) denotes a result obtained with the 
3PN-accurate EOB Hamiltonian using the LO (or OPN-accurate) spin-orbit terms, while E(3, 1) uses the 3PN- 
accurate EOB Hamiltonian with NLO (1PN- accurate) spin-orbit terms. Each panel in the Figure corresponds 
to a specific value of the dimensionless effective spin ciq. To guide the eye we use in all our figures a solid line 
to denote our "best" description, i.e. the 3PN-NLO EOB E(3, 1). Note that the addition of the NLO effects 
in the spin-orbit interaction has the clear effect of moderating the influence of the spins (especially for positive 
spins). While the binding energy curves using the LO spin-orbit effects tend to abruptly dive down towards 
very negative energies when the spins are large and positive, 16 the corresponding NLO curves have a much more 
moderate behaviour. 

Among the binding energy curves shown above, all the EOB ones (at least when the effective spin is not too 
large and positive), and some of the Taylor ones, admit a minimum for a certain value of the orbital frequency f2. 
This minimum corresponds to an inflection point in the corresponding (EOB or Taylor) Hamiltonian considered as a 



15 Indeed, in the physically most important case of parallel (rather than anti-parallel) spins, the spin-orbit coupling will be repulsive (like 
the effect of a positive 015), and will tend to reinforce the "bad" behaviour of the 3PN orbital Taylor Hamiltonian (i.e. the absence of 
any last stable orbit). 

16 As discussed in Section 3C of Ref. p[, this is due to the then repulsive character of the spin-orbit (and spin-spin) interaction. 
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FIG. 3: Energy e versus angular frequency Q along circular orbits for various values of the parameter fio, as predicted by the 
EOB Hamiltonian. We have assumed mi = mz, ai = «2, and 9o = n/2. As before E(n, s) refers to an EOB Hamiltonian, with 
nPN accuracy in the orbital terms, and an accuracy in the spin-orbit coupling equal to the LO one if s = 0, and the NLO one 
if s = 1. In all cases, we include the full LO spin-spin coupling. 



function of r. In other words, the minimum is the solution of the two equations 17 

dH Ica \ d 2 H 1 

cal 



dr 



0. 



dr 2 



0. 



(5.10) 



17 Note in passing that, in the EOB case, the two Eqs. I|5,10fl are equivalent to the two similar equations involving the effective Hamiltonian: 
dH eB /dr = 0, d 2 H cS /dr 2 =0. 
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TABLE I: LSO parameters for symmetric binary systems (with mi = mi and fii = 0,2 = ao) for the 3PN-NL0 EOB 
Hamiltonian E(3, 1). 



ao 


e 


Cl 


-1.00 


-0.01039 


0.04473 


— u. / 




u.uoioy 


-0.50 


-0.01270 


0.05989 


-0.25 


-0.01437 


0.07143 


0.00 


-0.01670 


0.08822 


0.25 


-0.02026 


0.11521 


0.50 


-0.02660 


0.16444 


0.75 


-0.03701 


0.23249 


1.00 


-0.03826 


0.22210 



The solutions of these two simultaneous equations correspond to what we shall call here the Last Stable (circular) Orbit 
(LSO). 18 Several methods have been considered in the literature 0, HH for using PN-expanded results to estimate 
the characteristics of the LSO. One of these methods consists in considering the minima in the Taylor expansion of 
the function e(Q). These minima (called "Innermost Circular Orbit" (ICO) in Refs. [H, [23j], where they were used to 
estimate the LSO of spinning binaries) differ from the minima in the Taylor energy curves considered in Fig. 1 above, 
which were based on using a Taylor-expanded Hamiltonian. The advantage of consistently working (as we do here) 
within a Hamiltonian formalism is that we are guaranteed that the minima in the corresponding energy curves, when 
they exist, do correspond to a Last Stable orbit (and an associated inflection point) for some well-defined underlying 
dynamics. By contrast the dynamical meaning (if any) of a minimum of the Taylor-expanded function e Taylor (f2) is 
unclear. Anyway, as we saw above that the 3PN-accurate Taylor-expanded orbital Hamiltonian does not admit any 
Last Stable Orbit, we have not plotted in Fig. 4 the Taylor-based predictions for spinning binaries because they do 
not seem to lead to reasonable results. 

Concerning the dynamical meaning of the LSO, let us recall that it had been analytically predicted in [|| (and con- 
firmed in recent numerical simulations that the transition between inspiral and plunge is smooth and progressive, 
so that the passage through the LSO is blurred. In spite of the inherent "fuzziness" in the definition of the LSO, it 
is still interesting to delineate its dynamical characteristics because they strongly influence some of the gross features 
of the GW signal emitted by coalescing binaries (such as the total emitted energy, and the frequency of maximal 
emission). 

Let us comment on the results of our study of the characteristics of LSO's: 

• In Fig. 4 we plot the LSO binding energy, predicted by the EOB approach, as a function of the dimensionlcss 
effective spin parameter do. We contrast LO spin-orbit versus NLO spin-orbit (1 versus 0). We use 3PN 
accuracy (for the orbital effects) in all cases, and always include the LO spin-spin interaction. The upper panel 
shows that the use of LO spin-orbit interactions leads to dramatically negative LSO binding energies when 
the spins become moderately large. [The middle panel is a close-up of the upper one, and focuses on spins 
So < +0.2.] We find that the 3PN-LO EOB Hamiltonian E(3, 0) admits an LSO only up to spins as large as: 
So < +0.9. However, as first found in [1], spin effects become dramatically (and suspiciously) large already 
when So > +0.5. By contrast, as we found above, the inclusion of NLO spin-orbit interactions has the effect 
of moderating the dynamical influence of high (positive) spins. The bottom panel focusses on our "best bet" 
3PN-NLO Hamiltonian E(3, 1). 

As mentioned above, Ref. [l5[ has considered, instead of the Taylor-Hamiltonian LSO, the minimum of the 
Taylor-expanded function e Taylor (17) (or "ICO"). For the two cases So = — 1 , (corresponding to their /t, = —1,0), 
they found, in the 3PN-NLO case, energy minima equal to e = E\co/ m = —0.0116, —0.0193 for corresponding 
orbital frequencies f2 = 77io;ico = 0.059,0.129. These numerical values should be compared with the numerical 
values we quote in Table I below. On the other hand, for the large and parallel spin case So = +1 Ref. [HI found 
that the Taylor-expanded function e Taylor (0) has no minimum. Finally, note that the qualitative shape of the 
curve giving the (EOB) LSO energy as a function of the spin parameters S is similar both to the corresponding 



As we recalled above, spinning binaries admit, in general, only spherical orbits, rather than circular ones. Reference []J studied the 
binding energies of the Last Stable Spherical Orbits (LSSO). Here, however, we restrict ourselves to the parallel spin, where it makes 
sense to study circular, equatorial orbits. 
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FIG. 4: Binding energy of the Last Stable (circular) Orbit (LSO) predicted by the EOB approach. We study the effect of 
including NLO spin-orbit terms by contrasting the LO and NLO predictions. We plot the dimensionless energy eLSO of the 
LSO versus ao- We have assumed mi = mi, ai = a,2, and 8o — n/2. For E(3, 0) a LSO exists up to So < +0.9. 



curve for a spinless test-particle in a Kerr background (see, e.g., Fig. 7 below), and to the curve giving the LSO 
energy of a spinning test particle in a Kerr background, as a function of the test spin (see Fig. 4 in Ref. [241]). 

To complement the information displayed in Figs. 1-4, we give in Table I the numerical values of the main LSO 
characteristics (binding energy and orbital frequency) for our "best bet" Hamiltonian, namely the 3PN-NLO EOB 
one E(3,l). 

In Figs. 5 and 6 we study the effective-spin-dependence of another LSO-related physical quantity of relevance for 
the dynamics of coalescing binaries: the total (orbital plus spin) angular momentum of the binary when it reaches 
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FIG. 5: The dimensionless total angular momentum Kerr parameter aj SO , Eq. (|5. 13|> . at the LSO, versus ao- We have assumed 
mi = ma, ai = as, and 0o = i~/2. The parameter aj SO is computed from Eq. (|5.13[l with j'lso = ^lso + Si + &2 = ^lso + 2oo- 
We compare the various EOB predictions obtained either by improving the accuracy of spin-orbit terms [E(3, 1) versus E(3, 0)], 
or by improving the accuracy of orbital terms [E(4, 1) versus E(3, 1)]. We use two representative values of the 4PN parameter 
as = +25 and as = +60. For comparison, we also include a fit to recent numerical estimates of the final Kerr parameter of the 
black hole resulting from the coalescence of the two constituent black holes. 



the LSO [i.e., at the end of the (approximately) adiabatic inspiral, just before the plunge], 

J = L + S! + S 2 . (5.11) 

In terms of rescaled dimensionless variables, this becomes 

; c - mi m? 

GM n mi mi 

where £ = c£. Actually, the most relevant quantity is the dimensionless Kerr parameter associated to the total LSO 
mass-energy and the total LSO angular momentum, i.e., the value at the LSO of the ratio 

aj = 0=" 9! (5. Id) 

G(iWc 2 ) (ffreal/(M C 2)) 2 

where j is the modulus of j. 

• In Fig. 5 we contrast the dependence of &j SO on the dimensionless effective spin parameter do for several EOB 
models: the two 3PN-accurate ones [E(3, 0) using LO-accurate spin-orbit, and E(3, 1) using NLO-accurate spin- 
orbit], and two illustrative @, "4PN-accurate" NLO-spin-orbit models E(4, 1) (using either — +25 or 
as = +60, as in Fig. 1). [Here, we are still considering fully symmetric systems with m\ — m% and a\ =02, 
so that do = a 1 — 0,2] Again we see the moderating influence of NLO corrections. The EOB-LO curve E(3, 0) 
exhibits a sudden drop down (pointed out in [l|) before rising up again (and disappearing at do = +0.9 when 
the LSO ceases to exist). By contrast, the NLO curve E(3, 1) exhibits a much more regular dependence on do, 
which is roughly linear over the entire range of values —1 < do < 1. The two illustrative E(4, 1) curves exhibit 
a "mixed" behaviour where a "drop" similar to the one featuring in the LO curve is still present, though it is 
moderated by NLO spin-orbit effects. This sensitivity to the inclusion of a 4PN contribution in A(u) is due to a 
delicate interplay between the modified shape of the basic spin-independent "radial potential" A{u, 05) and the 
use of a (1,4) Pade resummation of the "effective spin-dependent radial potential" A t (R), Eq. (|4.6a|) . Indeed, 
the additional contributions proportional to 05 and a 2 are both repulsive , and tend to compound their effect, 
which is to push the LSO toward closer, more bound orbits [l[. 

We have also indicated in Fig. 5 the final (i.e., after coalescence) dimensionless Kerr parameter of (symmetric) 
spinning binaries, as obtained in recent numerical simulations [26|, H3, Hg, 29]. For simplicity, we have shown the 
simple analytic fit proposed in 28]. The fact that the 3PN-NLO-accurate EOB LSO Kerr parameter [E(3, 1)] 
is systematically above the final Kerr parameter is in good agreement with the fact that, after reaching the 
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FIG. 6: The dimensionless total angular momentum Kerr parameter aj at the E(3,l) LSO versus an for various values of the 
parameter a\. Here we consider spin-dissymmetric systems with a\ ^ an (but still mi = mn and Go = 7r/2). The parameter 
&j SO is computed from Eq. (|5.13|l with j'lso = ^lso + 0,1 + 02. 
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FIG. 7: Comparison of the test-mass limit mn — /i — > (with fixed an; so that Sn/mn — * 0) for two Hamiltonians. We consider 
the specific "non relativistic" binding energy eLSO = e-LSo/v = (E-lso — Mc 2 )/(uc 2 ) at the LSO versus ao- The solid curve is 
the result of taking the test-mass limit of the EOB Hamiltonian, while the short-dashed curve is the result for a test particle 
moving in the Kerr spacetime. 



LSO, the system will still loose a significant amount of angular momentum 19 during the plunge and the merger- 
plus-ringdown. In the case of non-spinning binaries, it has been shown that, by using the EOB formalism up 
to the end of the process [i.e., by taking into account the losses of J and E during plunge, as well as during 
merger-plus-ringdown] , there was a good agreement (better than ~ 2%) between EOB and numerical relativity 
for the final spin parameter [l2T ] . We hope that the same type of agreement will hold also in the case of spinning 
binaries considered here. 

• In Fig. 6 we plot the LSO dimensionless Kerr parameter of Eq. (|5.13[) for spin- dissymmetric systems, namely 
at 7^ «2 (but with nil — 1^2), computed with the 3PN-NLO EOB Hamiltonian model E(3, 1). This plot 
illustrates that the LSO spin parameter is a smooth (and essentially linear) function of the two individual spins. 

• Finally, we compare in Fig. 7 the spinless test particle limit [i.e., mi — > 0, together with — S2/(m2c) — > 0, as 
appropriate to black holes for which a < 1] for two Hamiltonians: the 3PN-NLO EOB one E(3, 1), and the exact 



We use here the fact (found in numerical calculations, and implied by the analytical EOB approach), that, fractionally speaking, the 
angular momentum loss after the LSO is significantly higher than the corresponding energy loss. 
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one, as known from the geodesic action of a spinless test particle in the Kerr metric. For non-spinning systems 
the EOB Hamiltonian is constructed so as to reduce to the exact Schwarzschild-derived one in the test-particle 
limit. However, for spinning systems, we have chosen in Eq. (|4.6a|) to define the crucial metric coefficient A t (R) 
by Pade-resumming the sum of A(R; v) + a 2 / R 2 . This Pade-resummation is indeed useful for generally ensuring, 
for comparable mass systems, that At(R) have a simple zero at some "effective horizon" rjj. However, in the 
test-mass limit v — > 0, while the Taylor-approximant to A(R: v) + a 2 / R 2 would coincide with the exact Kerr 
answer, the Pade-resummed version of A(R\v) + a 2 /R 2 differs from it. We see, however, on Fig. 7 that the 
resulting difference has a very small effect on the LSO energy per unit (/x) mass, except when the dimensionless 
effective spin a is very close to +1. On the other hand, as we saw above when discussing Fig. 5, the issue of 
the Pade resummation of A t (R) becomes more subtle when one considers the comparable-mass case, together 
with the inclusion of a repulsive 4PN parameter a$ . 

VI. CONCLUSIONS 

The main conclusions of this work are: 

• We have prepared the ground for an accurate Effective One Body (EOB) description of the dynamics of binary 
systems made of spinning black holes by incorporating the recent computation of the next-to-leading order 
(NLO) spin-orbit interaction Hamiltonian [l3[ (see also Refs. [3 EH) into a previously developed extension of 
the EOB approach to spinning bodies [l[ . 

• We found that the inclusion of NLO spin-coupling terms has the quite significant result of moderating the effect 
of the LO spin-coupling, which would, by itself (as found in Ref. [lj), predict that the Last Stable (circular) 
Orbit (LSO) of parallely-fast-spinning black holes can reach very large binding energies of the order of 30% 
of the total rest-mass energy Mc 2 . By contrast, the inclusion of NLO spin-orbit terms predicts that the LSO 
of parallely-fast-spinning systems, though significantly more bound than that of non-spinning holes, can only 
reach binding energies of the order of 4% of the total rest-mass energy Mc 2 (see Fig. 4 above). This reduction 
in the influence of the spin-orbit coupling is due to the fact that the (effective) "gyro-gravitomagnetic ratios" 
are reduced by NLO effects from their LO values gj^° = 2, gj^? = | to the values (here considered along circular 
orbits) 

where x ~ GMj \Rc 2 ) ~ (GMfl/c 3 ) 2 ^ 3 . This reduction then reduces the repulsive effect of the spin-orbit 
coupling which is responsible for allowing the binary system to orbit on very close, and very bound, orbits (see 
discussion in Section 3C of Ref. [l|). 

• We studied the dependence of the dimensionless Kerr parameter of the binary system, dj = cJ / (G(H Ica i/ c 2 ) 2 ) , 
computed at the LSO, on the spins of the constituent black holes. Again the moderating effect of including 
NLO spin-orbit terms is very significant (compare the solid and the dashed 20 lines in Fig. 5). Thanks to this 
moderating effect the LSO Kerr parameter fij SO is found to have a monotonic, and roughly linear, dependence 
on the spin parameters of the individual black holes (see solid line in Fig. 5 and the various curves in Fig. 6). 
We also studied the effect of including the type of 4PN parameter a$ found useful in recent work 0, [1, HI [I3| 
for improving the agreement between EOB waveforms and numerical ones. 

• We leave to future work the analog of what was initiated for spinning systems in Ref. 5] , and recently completed 
for the case of non-spinning black holes in Ref. [12], i.e., a full dynamical study, within the EOB approach, of the 
Kerr parameter of the final black hole resulting from the merger of spinning black holes which takes into account 
the angular momentum losses that occur after the LSO, during the plunge, the merger, and the ringdown. Let 
us also note that Ref. [30j has recently proposed an approximate analytical approach (which is similar in spirit 



Compare also with Fig. 2 of Ref. fj| where the relevant LO result is the curve labelled "DJS" which reaches a maximum around 
a = ^ao — 0.31, in agreement with the (local) maximum in the dashed line of our Fig. 5 reached around ao ~ 0.36. 
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to the approximation used in Refs. [l|, d, [|| and above, namely that of considering the Kerr parameter of an 
effective test particle at, or after, the LSO) towards estimating the final spin of a binary black hole coalescence. 
The resulting prediction is, however, only in coarse agreement ~ 10% with numerical results. Note in this 
respect that, as displayed in Fig. 5, the "zeroth order" EOB result [corresponding to using the Kerr parameter 
for E(3, 1) at the LSO, without taking into account the later losses of angular momentum] is already in ~ 20% 
agreement with the fit to the numerical data [28|. The fact (displayed on Fig. 5) that the E(3,l) EOB LSO Kerr 
parameter is systematically above the final (after coalescence) Kerr parameter determined by recent numerical 
simulations [261 . 1271 . |28| . 1291 ] is in qualitative agreement with the fact that the system will loose a significant 
amount of angular momentum during the plunge and the merger-plus-ringdown. Note, however, the sensitivity 
of a^ so to a "4PN deformation" of the EOB Hamiltonian by the parameter a^. As said above, this sensitivity 
is due to the fact that the radial function A t (R)/R 2 combines the additional repulsive effects of both a positive 
4PN contribution +a$v{GM / '(c 2 R)) 5 and a positive spin-dependent contribution +a 2 /R 2 . We leave to future 
work an exploration of this issue, which might need the use of a different Pade resummation than the (1,4) one 
used in (|4.6a[) . 

It remains to be seen whether the EOB/Numerical Relativity comparison for the final Kerr parameter of spinning 
systems will be as good as it was found to be for the non-spinning case [Hj], i.e., at the 2% level. If this is the 
case, it will establish the physical relevance of the improved EOB Hamiltonian constructed in the present paper. 

• Let us finally note that there is some flexibility in the improved spin-dependent EOB Hamiltonian proposed 
above (besides the flexibility in the choice of the Pade resummation mentioned above). On the one hand, the 
choice p. 141) for the gauge parameters a(y) and b(v) might be replaced by other choices. On the other hand, 
the choice (14. lip for the effective spin vector might also be replaced by other ones. In particular, it might be 
interesting to consider the alternative definition 

A/ca ncw = S e ffnew = 

= ^s ffnow (S + S*). (6.2) 

This definition coincides with the one used above at LO in spin- orbit effects (because gf ncw = 2 + 0(v/c 2 )), 
and allows one to use a simplified supplementary spin-orbit contribution, built with 

a^^^gf.-gf)^, (6.3) 

instead of (|4.15|) . It might be interesting to explore which of these possible definitions exhibits the best agreement 
with current numerical results. 
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